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LIMITED-ANGLE FREQUENCY- DISTANCE RESOLUTION RECOVERY IN 

NUCLEAR MEDICINE IMAGING 

Background of the Invention 

The present invention relates to the diagnostic 
imaging art. It finds particular application in 

conjunction with nuclear single photon emission computed 
5 tomography (SPECT) medical imaging and will be described 
with particular reference thereto. However, the invention 
will also find application in conjunction with other types 
of non- invasive diagnostic imaging. 

Heretofore, diagnostic images have been 

10 generated from single and multiple-head nuclear cameras. 
Typically, a patient positioned in an examination region 
is injected with a radiopharmaceutical. Heads of the 
nuclear camera are positioned closely adjacent to the 
patient to monitor the radiopharmaceutical. Typically, 

15 the heads are stepped in increments of a few degrees 
around the patient until 360** of data have been acquired. 
That is, projection data along directions spanning 360" 
are collected. With multiple-head systems, projections 
along each direction need only be collected with one of 

20 the heads and be assembled into a complete data set. 

Each detector head carries a collimator which 
defines a path along which it can receive radiation. 
However, due to the finite length and dimensions of the 
collimator, each incremental area of the detector head 

25 actually views an expanding cone. Thus, with increasing 
depth into the patient away from the detector head, the 
region from which a sensed radiation event originated 
expands. This creates depth-dependent blurring and 



- 2 - PKR 2 0 655 

uncertainty in the resultant image data. This error is a 
non- stationary convolution which is difficult to 
deconvolve. However, when the patient is viewed over a 
full 360", the angular data sets are periodic in 27r 
5 radians. By transf orming^ the full data sets into the 
frequency domain with a Fourier transform fit to the 
sampling intervals, the non- stationary deconvolution is 
reduced to a stationary deconvolution problem, 
particularly for high frequencies. 

10 Although these prior art resolution recovery 

techniques work well on full data sets, cardiac imaging is 
typically done using an incomplete data set. More 
specifically, in a three-head camera system where only two 
heads collect the emission data while the third head is 

15 used to collect transmission data, the gantry is rotated 
by about 102" to generate the equivalent of only about 
204' of emission data. The prior art resolution recovery 
techniques do not work on partial data sets whose data is 
not periodic in 27r radians. 

2 0 The present invention contemplates a new and 

improved method and apparatus which overcomes the 
above -referenced problems and others. 

Summary of the Invention 

25 In accordance with one aspect of the present 

invention, a method of diagnostic imaging is disclosed. 
A plurality of projection data sets are collected at each 
of a plurality of angles around a subject. The projection 
images are collected over less than 360". A resolution 

3 0 recovery process is performed on the projection data sets. 

The resolution recovered projection data sets are 
reconstructed into an image representation. 

In accordance with another aspect of the present 
inyention, a method of diagnostic imaging is disclosed. 
3 5 A gantry moves a detector head in a continuous angular 
orbit about a subject in an examination region. Data is 
collected during the continuous orbit and sorted into a 
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plurality of projection data sets corresponding to each of 
a plurality of angular increments around a subject. A 
resolution recovery process is performed on the projection 
data sets. The resolution recovered projection data sets 
5 are reconstructed into an image representation. 

In accordance with yet another aspect of the 
present invention, a diagnostic imaging apparatus is 
disclosed. At least one detector head detects incident 
radiation. A collimator mounted to the detector head 

10 limits trajectories along which radiation is receivable by 
the head. A movable gantry moves the detector head around 
a subject in an examination region. A data acquisition 
system acquires projection data sets from the detector 
head at angular increments spanning less than 360". A 

15 zero-filling processor generates zero-filled data sets 
between the actually collected projection data sets, to 
create 360* of data sets. A smoothing processor smooths 
interfaces between the actually collected and zero-filled 
data sets. A resolution recovery processor operates on 

20 the smoothed data sets. A reconstruction processor 
reconstructs the resolution recovered data sets into a 
three-dimensional image representation. An image memory 
stores the three-dimensional image representation. 

One advantage of the present invention is that 

25 it accurately restores limited-angle data sets. 

Another advantage of the present invention is 
that it restores continuously scanned data sets. 

Another advantage of the present invention is 
that it processes image data in a clinically feasible 

30 time. 

Still further advantages and benefits of the 
present invention will become apparent to those of 
ordinary skill in the art upon reading and understanding 
the following detailed description of the preferred 
3 5 embodiment s . 
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Brief Description of the Drawings 

The invention may take form in various 
components and arrangements of components, and in various 
steps and arrangements of steps. The drawings are only 
5 for purposes of illustrating preferred embodiments and are 
not to be construed as limiting the invention. 

FIGURE 1 is a diagrammatic illustration of a 
nuclear medicine imaging system in accordance with the 
present invention; 
10 FIGURE 2 is a diagrammatic illustration of the 

formation of the modified data set which is periodic in 2tt 
radians, based on the non-periodic limited- angle data set; 

FIGURE 3 shows the coordinate system used in 
data acquisition and image reconstruction. 

15 Detailed Description of the Preferred Embodiments 

With reference to FIGURE 1, a nuclear camera 
system 10 includes a plurality of detector heads 12, in 
the preferred embodiment three detector heads 12^, 123, and 
123. Two of the three heads are typically used to obtain 

20 emission data, while the third head is used to obtain 
transmission data. Of course, other numbers of detector 
heads can also be utilized. Each of the detector heads 
includes a collimator 14^, 142, and 143. In the preferred 
embodiment, the collimators collimate incoming radiation 

25 from a subject 16 to parallel rays. However, because the 
collimators have finite size, each collimator permits rays 
which lie along a corresponding cone to pass to the 
detector head. The cone expands with depth into the 
patient from the detector head. 

30 Each of the detector heads includes a 

reconstruction system 20-l, 20^, and 2O3 which determines the 
coordinates on a face of the detector head in the 
longitudinal or z-direction of the patient and the 
transverse direction across the detector head. With the 

35 detector head in a single orientation, scintillation 
events are collected for a preselected period of time to 
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generate a projection image representation. After the 
preselected data acquisition duration, a rotating gantry 
22 rotates all three detector heads concurrently a short 
angular distance, e.g., 3*. In the new location, each of 
5 the detector heads collects another projection image. An 
angular orientation monitor 24 determines the angular 
orientation of each of the heads at each angular data 
collection position. 

A data acquisition system 30 receives each of 

10 the projection images and an indication of the angle along 
which it had been taken. This data is stored in a 
three-dimensional memory 32 in longitudinal (z) , lateral 
(s) , and angular (0) coordinates. 

With continuing reference to FIGURE 1 and 

15 further reference to FIGURE 2, the rotating gantry 22 
rotates the detector heads over 102* for cardiac imaging. 
Rotation over 102' generates the equivalent of about 204* 
of emission data 110 collected by two camera heads. Each 
data set is collected at equal angular increments, 

20 e.g., 3**. FIGURE 2 shows an incomplete data set with no 
data collected between 204** and 360* in the angular 
dimension. Simultaneously, the third camera collects 
transmission data over the 102* gantry rotation. 

With continuing reference to FIGURE 2 and 

25 reference again to FIGURE 1, a resolution recovery system 
42 includes an angular displacement dimension resolution 
recovery system and optionally includes resolution 
enhancement sub- system 46 for the axial -dimension. The 
resolution recovery system in the angular direction (0) 

30 has a zero-filling processor 50 which creates zero 
magnitude projection data sets 112 at each of the 3* 
intervals between 204° and 360*. In this manner, a 
function which is periodic in 27r radians is created. In 
order to prevent Gibbs' ringing, the sudden 

3 5 discontinuities at each end of the actually collected data 
110 between the actually collected and zero-filled data 
112 are smoothed by a smoothing processor 56. In the 
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preferred embodiment, the magnitude of the end points are 
each cut in half producing modified data points 114 and 
116. Optionally, other smoothing functions which span 
several points at each end are also contemplated. 
5 A Fourier transform processor 60 transforms the 

data into the frequency domain. The Fourier transform is 
selected to match all of the sampling points including the 
zero- filled points in the angular dimension. A stationary 
deconvolution processor 62 operates on the 

10 frequency-spaced data. In frequency space, the 

deconvolution problem reduces to a stationary 
deconvolution problem, particularly for high frequencies. 
Because high frequencies correspond to fine detail, it is 
the high frequencies which correspond to the resolution to 

15 be optimized. The additional deconvolution 46 may 
optionally also be performed to improve image resolution 
in the z-direction. After the data has been deconvolved, 
an inverse Fourier transform processor 64 transforms the 
data from the frequency domain back into real space. 

20 Optionally, a three-dimensional memory 66 stores the 
resolution recovered sets of projection data. A 
reconstruction processor 68 reconstructs the projection 
data sets, using filtered back-projection, iterative 
reconstruction, or other techniques which are well-known 

25 in the art, and stores the resultant image in a 
three-dimensional image memory 70. A video processor 72 
withdraws selected portions of the reconstructed image and 
converts it to .appropriate format for display on a 
human- readable monitor 74, such as a video monitor, CCD 

30 display, active matrix, or the like. The video processor 
may withdraw selected slices', three-dimensional 
renderings, projection views, or the like. 

In one alternate embodiment, the movable gantry 
22 rotates the detector heads continuously. Although the 

35 detector heads rotate continuously, the data acquisition 
system 30 bins the collected data into regular angular 
intervals, e.g., 3'. Because the data is collected over 
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3' of rotation, there is an additional blurring component. 
As described in greater detail below, the stationary 
deconvolution processor 62 deconvolves the Fourier space 
data with respect to the blurring caused by the continuous 
5 motion and treats data collected over a few degrees as if 
it were all collected at precisely the same angle. 

With reference to FIGURE 3, and looking now to 
details of the resolution recovery system, the working 
principles are explained in two-dimensions and circular 

10 orbits, but the generalization to three-dimensions is 
straightforward and generalization to non-circular orbits 
is known to those conversed in the art. The third 
dimension, z, is perpendicular to the plane defined by the 
x-y axes. The scanner collects data in coordinates (s,0) 

15 which are the sinogram coordinates of the Radon transform 
coordinates. The unblurred or undegraded Radon transform 
of the object o is the line integral along t, the axis 
perpendicular to s. In a nuclear medicine tomographic 
device, the object o is blurred by a point response 

20 function g. This process is modeled by a convolution in 
s. The amount of blurring depends on the depth t. Hence, 
the blurring function is a non-stationary convolution of 
o with g. The result is the blurred Radon transform p. 
For limited-angle tomography, this operation of the 

25 scanner is represented by: 

oo oo 

p(s,<t)) = f j'g{s-s',t)o{s'e + te'-)ds'dt (i) , 

where d and 6^ are two-dimensional unit vectors aligned 
with the (s,t) axes: 

9 : [cos((), sinct)] and 9"^ : [ -sine}), coscj)] (2). 

A point source located at x in the object will have the 
point -source projection: 
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p^ = j^q {s-s\ t) \J\5 (s^-x-e, t-x-e-") ds^dt = g ( s-x- 0, x-e-") (3) , 


where |j| is the Jacobian of the change of coordinates 
from (x,y) to (s,t) and is equal to one because it is a 
rotation. The two-dimensional Fourier transform of 
Equation (3) gives: 


2n 


Pg(co,n) = Je~^'^*Je"^"^g(s-x-0,x-0-')dsd(t> 

0 

2n 

= |e'^"*e"^"^"®g(co,x-e^)dc{) ^ g 


CO, x-0-^ = 

CO 


(4) 


5 The </)-integral is evaluated using the principle of 
stationary phase, an approximation process familiar to 
those knowledgeable in the art, which asymptotically 
converges to the exact value of the p^ integral at high 
frequencies. In Equation (4), it will be noted that the 

10 solution does not depend on the location of the point 
source. Hence, it can be applied to any point source or 
a weighted collection of point sources. 

Furthermore, if the scanning arc is limited to 
something less than 360" or 27r radians > the basic result 

15 of stationary phase, that the distance t is given by: 

t = x-0^ = — (5) 


is unchanged, because every distance is represented in 
each planar projection image. This enables the limits of 
integration in Equation (4) to be replaced with an arc or 
a series of disjoint arcs and Equation (5) still holds to 
2 0 good approxima t ion . 
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In operation, the digitized three-dimensional 
data set p(Sj,Zj^,0^) is collected. The digitization assumes 
equal increments in each coordinate sample space. For 
example, the angular increment m is taken at intervals 
5 0^=27rm/N radians, where N is an integer. 

The depth dependent , point response function 
g(s,z,t) is obtained. When the scan is over a limited 
angle arc described by the interval [0u/0v]/ ^ digital 
input data set is formed according to the following 
10 process: 

0, for 0:^ W < min {u,v) 
— p(s^.,Zj^,ct)^) , for W = u 
p.^(s^., z^,ct)p^) = p(s^,z^,(t)p^)for u < iV< V (6). 
-|p(s^,z^,c|>^) , jfor IV = V 
0 , for max ( u, v) < W < N 


The endpoint projection frames are multiplied by one-half 
to prevent excessive ringing or Gibbs' phenomena. Other 
smoothing functions for smoothing the sudden data 
discontinuity are also contemplated. The Fourier 

15 transform with respect to angle is selected to have 
dimension N, no more and no less, to preserve the cyclic 
nature of the data. The three-dimensional Fourier 
transform applied by the Fourier transform processor 60 of 
the modified projection data p^^ is denoted by: 

Pin ("s'^z'^) (7) . 


20 


The filter function: 


CO , CO , 


n 

CO 


(8) 


is defined. The filter function is used in a regularized 
inverse filter such as a noise reduction filter or a Metz 
filter to perform the stationary deconvolution 62. The 
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inverse Fourier transform processor 64 takes the inverse 
Fourier transform of the filtered image to obtain 
projection data that possesses improved resolution. 

In the continuous scanning embodiment, the 
5 camera moves according to a preprogrammed orbit. In the 
preferred embodiment, the three heads move simultaneously 
in a complex patient dependent orbit. Continuous motion 
is advantageous for reducing motion complexity. However, 
continuous scanning degrades resolution. Preferably, the 
10 stationary deconvolution module includes a component to 
compensate for this degraded resolution. 

Continuous scan projection data is related to 
the step-and- shoot projection data p by: 

p (s,z,(J)) = f p{Sr z,(p^)d<p^ (9). 

^ Ad) J 

^ -A(D/2 

The Fourier transform relationship is: 

Pc(*^s'"z'^) = p(co , CO ,n) \ J ^ , - (10). 

^ ^ ^ s z nAd)/2 2 2 


15 These relationships hold true for non-circular and 
limited-angle orbits as well. The filter function for 
continuous sampling g^, is modified to include the weighted 
sine function from Equation (10) : 


n 

CO 


^ sin(nA(l)/2) J 


nA(|)/2 


n 

CO 


(11) . 


In accordance with another alternate embodiment, 
20 in the step- and- shoot imaging mode, the actually collected 
data integrals may be disjoint. It is not necessary for 
all of the angular views to be adjacent as in the example 
above . 

The invention has been described with reference 
25 to the preferred embodiments. Obviously, modifications 
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and alterations will occur to others upon reading and 
understanding the preceding detailed description. It is 
intended that the invention be construed as including all 
such modifications and alterations insofar as they come 
5 within the scope of the appended claims or the equivalents 
thereof . 


